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Background 

The accuracy of ligand-protein docking may be affected 
by the presence of water molecules on the surface of 
proteins. Water can form complex bridging networks 
and can play a critical role in dictating the binding 
mode of ligands. A recent analysis of high-resolution 
crystal structures of ligand-protein complexes revealed 
that 85% of the complexes had one or more water mole- 
cules bridging the interaction between ligand and pro- 
tein. For predicting the binding modes and energies of 
protein-ligand interactions, molecular docking methods 
are commonly used. In order to obtain an accurate 
complex geometry and binding energy estimation, an 
appropriate method for calculating partial charges is 
essential. AutoDockTools software, widely used as inter- 
face for preparing input files, utilizes the either Gastei- 
ger or Kollman partial charge calculation method for 
both protein and ligand charge calculations. However, it 
has already been reported that more accurate partial 
charge calculation and as a consequence, more accurate 
docking can be achieved by using quantum chemical 
methods. In common practice so far, the quantum che- 
mical partial charges were applied to the ligands for 
docking calculations. The newly developed Mozyme 
function of MOPAC2009 allows fast partial charge cal- 
culation of proteins by quantum mechanical semi- 
empirical methods. Thus, in the current study, we use 
the semi-empirical quantum-mechanical partial charge 
calculations to investigate the interaction energies and 
polarization effects of the various components of the 
binding pocket on a set of Cathepsin B protein. 



Results 

The docking accuracy was computed by using the origi- 
nal AutoDock scoring function with the set of 19 pro- 
tein ligand complexes using Gasteiger, AMI and PM3 
partial charge calculation methods. This helped us to 
compare the effect of the partial charge calculation 
method on docking accuracy. It was seen that the dock- 
ing accuracy in regard to complex geometry significantly 
increased when partial charges of the ligands and pro- 
teins were calculated with the semi-empirical PM3 
method. Our results demonstrate that (i) the energetic 
of the key water molecule are more favorable for the 
binding site in the Cathepsin B protein (ii) Water brid- 
ging and triangle formation were seen between the key 
amino acid residue and the ligand (iii) The internal 
energy is significant factor for the binding modes of var- 
ious ligands. It was also observed a statistically signifi- 
cant overall increase in accuracy when water molecules 
are included during docking simulations. Out of the 19 
complexes analyzed in the course of our study, the geo- 
metry of 17 complexes were accurately calculated using 
PM3 partial charges, while the use of Gasteiger charges 
resulted in only 8 accurate geometries. 

Conclusion 

Our findings indicate that the inclusion of water 
molecules in ligand-protein docking results in signifi- 
cant increases in docking accuracy when the use of 
quantum chemical partial charge assignment on both 
ligand and proteins for predicting the docking 
simulations. 
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(a)EPO with AM1 Charges 



(b)EPOwith AM 1 Charges 

HUMO surface 





(c) EPO with PM3 Charges 



Figure 1 Charge and force field applied to inhibitors for better docking score. 



(d) EPO with PM3 Charges 
HUMO surface 
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